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Our model of shaken sand, presented in earlier work, has been extended to include a 
more realistic 'glassy' state, i.e., when the sandbox is shaken at very low intensities of 
vibration. We revisit some of our earlier results, and compare them with our new results 
on the revised model. Our analysis of the glassy dynamics in our model shows that a 
variety of ground states is obtained; these fall in two categories, which we argue are 
representative of regular and irregular packings. 



1. Introduction 

The test of a good lattice model of a complex system is whether it succeeds in 
capturing the essential physics of a real system in its bid to reduce its technical 



complexity. Areas as diverse as traffic flow , plate tectonics [ 12 , Ij and granular 
flow are examples where lattice models have been used rather successfully 
despite their apparent simplicity, to describe at least a few of the salient features 
of some genuinely complex systems. 

In this spirit, we present two versions of a model of shaken sand in the following; 
both models exhibit behaviour that is representative of shaken sand between the 
fluidised and the glassy regimes. While the second model includes rather complex 
interactions between the 'grains' in the frozen state, in contrast to the first, the latter 
is nevertheless surprisingly successful in replicating at least some of the qualitative 
features associated with the glassy regime. One of the aims of this contribution 
is then to identify some of the essential features that are needed for a (discrete) 
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minimal model of the system in question. 

The earlier model ('old model') is the generalisation of a cellular-automaton 
(CA) model [p|j2^ of an avalanching sandpile. This model shows both fast and slow 
dynamics in the appropriate regimes: in particular, it reduces to an exactly solvable 
model of noninteracting grains in the frozen ('glassy') regime, and provides one with 
a toy model for ageing in vibrated sand |^J2^ . 

Next, we present a more realistic version of the model ('new model') which 
is the topic of ongoing research. While identical in the fluidised regime, the latter 
model is less of a toy model for the glassy state, in that the grains are no longer non- 
interacting, but are coupled to each other based on their orientations. Our analysis 
shows that a rich variety of ground states is obtained, which we analyse in terms 
of a particular parameter that has an interpretation in terms of the irregularity of 
the grains. 

2. The Model 

We consider a rectangular lattice of height H and width W with N < HW grains 
located at the lattice points, shaken with vibration intensity F. Each grain is a 
rectangle with sides 1 and a < 1, respectively. Consider a grain in row i and 
column j whose height at any given time is given by hij = riij^ + auijj^, where 
Uij- is the number of vertical grains and Uij^ is the number of horizontal grains 
below (j, j). 

The dynamics of the model are described by the following rules: 

(i) If lattice sites (i + 1, j — 1), [i + or (i + l,j + 1) are empty, grain 
(i,j) moves there with a probability exp(— l/F), in units such that the 
acceleration due to gravity, the mass of a grain, and the height of a lattice 
cell all equal unity. 

(ii) If the lattice site {i — 1, j) below the grain is empty, it will fall down. 

(iii) If lattice sites {i — l,i ± 1) are empty, the grain at height hij will fall to 
either lower neighbour, provided the height difference hij — /ii_i,j±i > 2. 

(iv) The grain flips from horizontal to vertical with probability exp(— m.y (Ai7-|- 
A/i)/r), where rriij is the mass of the pile (consisting of grains of unit 
mass) above grain For a rectangular grain, A_ff = 1 — a is the height 
difference between the initial horizontal and the final vertical state of the 
grain. Similarly, the activation energy for a flip reads A/i = 5—1, where 
b = vl + c? is the diagonal length of a grain. 

(v) The grain flips from vertical to horizontal with probability exp{—mijAh/T). 

These rules yield a rich variety of dynamical behaviour, depicted in Fig. |l|, which 
shows the time evolution of the sandbox under the effect of different vibrational in- 
tensities. All the boxes are started in the same initial state; the observed behaviours 
correspond to the glassy regime (top row) and the fluidised regime (bottom row). 
In the intermediate regime separating the two p2|, individual-particle relaxations 
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Fig. 1. Snapshots of a sandbox with Vy = 30, H = 100, n = 2500 grains with a = 0.7, A?t = 0.05, 
and shaking intensities F = 0.1 (top row: glassy regime) and 0.8 (bottom row: fiuidised regime), 
for times t = 1 (before tapping), t = 2 (each grain on average touched once by the MC simulation), 
i = 5, 10, 30, and 100. 



are the mechanism for each particle to find local stability, fn related work 

this was associated with a threshold, called the single particle relaxation threshold 

(SPRT). 

In line with recent investigations of compaction ||l|,|,|,0|l|,|l|,|o|,|7l|^,|l| , we 
have examined the behaviour of the packing fraction of our model, as a function of 
the vibration intensity F. Let N~ and iV+ be the numbers of vertical and horizontal 
grains in the box. The packing fraction cj) is: 

, N+-aN- 

— \ , 2.1) 

which we use as an order parameter reflective of the behaviour of the compactiv- 
ity |l9| , po| . The vertical orientation of a grain thus wastes space proportional to 1— a, 
relative to the horizontal one. 

The advantage of our current model is that at least conceptually it can be 
extended to non-rectangular grain shapes. In this sense, we can and will consider a, 
AT?, and A/i as phenomenological parameters in the following. 
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3. A Spin-Model for the Ordered Regime 

The frozen regime is characterised by an absence of holes within the sandbox, and 
negligible surface roughness. Here, the earlier model reduces to an exactly 
solvable model of W independent columns of H noninteracting grain orientations 
o'nit) = ±1, with a — +1 denoting a horizontal grain, and a — —1 denoting a 
vertical grain. The orientation of the grain at depth n, measured from the top of 
the system, evolves according to a continuous-time Markov dynamics, with rates 

rw;(- ^+) =exp(-nA/i/r), 

\ wi+ ^ -) - exp(-n(Ai? + Ah) /F), ^ ' 

as rriij = ?i = H+l—i. The parameters AH and Ah correspond to two characteristic 
lengths of the model, 

^eq = r/AH, ^dyn = T / Ah. (3.2) 

The equilibrium length is the typical depth below which all grains are frozen 
into their horizontal ground-state orientation, while the dynamical length ^dyn is 
characteristic of the divergence of the relaxation time with depth, r ~ exp(n/^dyn)- 
As a consequence, any perturbation propagates only logarithmically slowly down 
the system, over an ordering length A{t) w ^dyn Int. 

The detailed analysis of this model in earlier work, despite its 'toy model' nature, 
was remarkably successful in reproducing certain features of the glassy state. In this 
work, we present an improved version where the spins are no longer noninteracting, 
in the frozen regime. They are in fact coupled in a rather complex way; these 
more realistic interactions lead to a rich variety of ground states depending on the 
analogue of the aspect ratio a, which we argue is representative of the shape of the 
grains. We present this model in the next section. 



4. A Generalisation 

The generalisation of our earlier model involves the insertion of eq. (2.1) into the 
transition rates of the system. We thus require that, for a given value of a, the 
transitions are such that the packing fraction of the system is locally minimised. 
We now allow a to take arbitrary values (while always remaining positive): I — a can 
then be visualised as the size of the 'void' associated with the 'wrong' orientation 
of the grain. For convenience, we write these rules in a more general form below. 

Thus, the dynamical rules (iv) and (v) of the rectangular grain model may be 
regarded as a special case of a more general model with transition rates 

r w{ >+)= exp(-(Ay + r7ij)/r): 

l«;(+^-)=exp(-(A,,-r7,,)/r), 

where 

ri.j = Am,j+ + BrUij^, = Cmi-j+ -f Dm,j^, (4.2) 
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are, respectively, the ordering field and the activation energy felt by grain {ij). In 
these expressions, mij± is the number of horizontal and vertical grains above the 
grain respectively. Note that grains are only counted from to the first 

void above level j. Thus my = my+ + rriij-. 
The earlier model is recovered by setting 

A = B = AH/2, C ^ D = AH/2 + Ah. (4.3) 

In the general situation where the equalities A = B, C = D are not obeyed, 



the rates (4.1) depend on the orientations of all the grains above the grain under 
consideration. Our new model is therefore a fully directed model of interacting 
grains, where causality acts both in time and in space, as the orientation of a given 
grain only influences the grains below it and at later times. 

The key parameter which governs the statics and dynamics of the model at low 
shaking intensity turns out to be the dimensionless ratio e ~ A/B. Consider the 
ordered regime, where there are no holes, in the zero-temperature limit (F — + 0). 
In this regime, the steady-state values of the grain orientations are given by the 
deterministic, recursive equation 

o-„ = sign(£n+(n) + n_(n)), (4.4) 

where n±{n) is the number of horizontal and vertical grains at depths 0, . . . , n — 1, 
so that n+(n) + n^(n) ~ n = 1,2, . . . We assume ctq = +1. 



As long as £ > 0, (4.4) leads to the trivial ground state where all the grains are 
horizontal ((t„ = -|-1), generalising thus the case of the earlier model (e = 1). 

In the frustrated regime (e < 0), ground states have a richer structure. They 
contain non-trivial fractions of horizontal and vertical grains, 

/+ = l/(l-e), /_ = -e/(l-e). (4.5) 

The grain orientations are distributed in a way which depends, rather unexpectedly, 
on whether s is rational or not. 

• Rational case: if e = ^p/q is a negative rational number (in irreducible 



form), (4^) determines the grain orientation cr„ whenever the depth n is 
not a multiple of r = p+q. The orientations of the latter grains are left free. 
We thus obtain an extensively degenerate set of ground states, each of them 
being a random sequence of two types of r-mers, i.e., clusters of r grains. 
The associated configurational entropy per grain reads S = (ln2)/r. The 
simplest example is e = — 1, hence r = 2 and /_ = 1/2, where the clusters 
are H — and — + , so that the ground states are all the dimerised grain 
configurations. Two examples consist of trimers (r — 3), namely e — —1/2, 

with /_ = 1/3 and clusters H \- and — h +, and e — —2, with /_ = 2/3 

and clusters H and — I — . 



Irrational case: if e is a negative irrational number, (4.4) determines all 



the grain orientations, so that the model admits a unique, non-degenerate 
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ground state, where the grain orientations are distributed in a quasiperi- 
odic fashion. The rule (4.4) is indeed equivalent to the cut-and-project algo- 
rithm used to build quasiperiodic binary chains, which are one-dimensional 
analogues of perfect quasicrystals ||l5|,|l7 18 2l| , p2| , p3|j2^ . For instance, for 
-s = (V5- l)/2 « 0.618033 (the inverse golden mean) , the unique ground 
state of the model is given by the well-known Fibonacci sequence: 



+ - 



5. Numerical Results 

In Fig. H we show the behaviour of the packing fraction as a function of F, the 
shaking intensity, for both models described above for a square box of side 120, 
containing 10000 grains of unit mass, each of which has an aspect ratio of 0.7. 

We note that the overall behaviour of the two models is rather similar, although 
the complexity of the second makes it far more prone to fluctuations even after the 
steady state has been reached. In each case, we observe 

• a fluidised region: for F 3> 1, we observe an initial increase (caused by a 
non- equilibrium and transient 'ordering' of grains in the boundary layer) of 
the packing fraction that quickly relaxes to the equilibrium values 0oo in 
each case. This over-shooting effect in Fig. || increases with F, since grains 
ever deeper in the sandbox can now overcome their activation energy to 
relax to the horizontal. 

• an intermediate region (for F w 1), where the packing fraction remains ap- 
proximately constant in the bulk, while the surface equilibrates via the fast 
dynamics of single-particle relaxation. The specific <^oo at which this occurs, 
is the single-particle relaxation threshold density observed in Ref. ^ ; non- 
equilibrium, non-ergodic, fast dynamics allows single particles locally to 
find their equilibrium configurations at this density. Analogous effects have 
been observed in recent experiments on colloids fssj] , where the correlated 
dynamics of fast particles was seen to be responsible for most relaxational 
behaviour before the onset of the glass transition. 

• a frozen region (for F <C 1), where the slow dynamics of the system results 
in a logarithmic growth of packing fraction with time: 

0^ = 6(F)lnt-ha, (5.1) 

where 6(F) increases with F, in good agreement with experiment ]30| , ^l]| . 
The slow dynamics has been identified § with a cascade process, where 
the free volume released by the relaxation of one or more grains allows 
for the ongoing relaxation of other grains in an extended neighbourhood. 
It includes the phenomenon of bridge collapse, which, for low vibration 
intensities, has been seen to be a major mechanism of compaction [p[|^,p7||. 
As F decreases, the corresponding (f>ac increases asymptotically towards the 
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Fig. 2. Behaviour of packing fraction for both models described in text, with an aspect ratio of 
0.7 in each case. For the old model (left), A = B = 0.35, C = D = 0.4, while for the new model 
(right) A = —0.2, B = C = D = 0.4. The shaking intensities are, for bottom to top, F = 0.3 (old 
model only), 0.5, 0.8, 1, 1.5, 2, 5, and 10. 



jamming limit ^jam, identified with a dynamical phase transition in related 
work 

While we have presented in earlier work a full analysis of two-time correla- 
tion functions for the old model, work is currently in progress to investigate this in 
the rather more complex new model. 

We finally investigate the analogue of 'annealed cooling', where F is increased 
and decreased cyclically, and the response of the packing fraction observed . 
The results obtained here are similar to those Q seen using more realistic models 
of shaken spheres, but the simplicity of the present lattice-based models allows for 
a greater transparency. 

Starting with the sand in a fluidised state, as in the experiment ||30| , ^ , we 
submit the sandbox to taps at a given intensity F for a time itap and increase the 
intensity in steps of (5F; at a certain point, the cycle is reversed, to go from higher to 
lower intensities. The entire process is then iterated twice. Fig. ^ shows the resulting 
behaviour of the volume fraction (/> as a function of F, where an 'irreversible' branch 
and a 'reversible' branch of the compaction curve are seen, which meet at the 
'irreversibility point' F* ||30| , |3l| . The left- and right-hand side of Fig. [| correspond 
respectively to high and low values of the 'ramp rate' (5F/ttap p0|j31[ , while the 
upper and lower panels correspond respectively to the old and new versions of our 
model. As the ramp rate is lowered, we note that: 

• the width of the hysteresis loop in the so-called reversible branch decreases, 
in both cases. The 'reversible' branch is thus not reversible at all; more 
realistic simulations of shaken spheres H,^,^ confirm the first-order, irre- 
versible nature of the transition, which allows the density to attain values 
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Fig. 3. Hysteresis curves. Top row: old model. Bottom row: new model. Left: 5r = 0.1, ttap = 
2000 time units. Right: ST = 0.001, ttap = 10^ time units. Note the approach of the irreversibility 
point r* to the 'shoulder' Fjam) the ramp rate ST /ttap is lowered. The packing fraction tends to 
its close-packing limit in the limit of low intensities for the old model, while it asymptotes towards 
the jamming limit for the new model (see text). 



that are substantially higher than random close packing, and quite close to 
the crystalline limit |^ . Precisely such a transition has also recently been 
observed experimentally in the compaction of rods . 
In both panels, the 'irreversibility point' T* approaches Fjam (the shaking 
intensity at which the jamming limit ^jam is approached), in agreement 
with results on other discrete models [Q. However, in the upper panel, 
the packing fraction at low intensities tends towards close-packing (the so- 
called dynamical transition referred to in |8j); this is at odds with the results 
of real experiment, which models with a greater degree of complexity 
have been able to replicate. In the lower panel, which corresponds to our 
new model, we see tentative indications of an improvement in this respect 
vis-a-vis the old model; the packing fraction here asymptotes towards the 
jamming limit, rather than rising indefinitely towards the close-packing 
limit 1,0. 



6. Discussion 

We have presented in the above two models of shaken sandboxes; while their design 
was such that they would show identical behaviour in regimes where there were a 
finite density of voids, the modelling of the densely packed regime was completely 
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distinct in each case. In the first case, the model reduced to a model of noninteract- 
ing spins, while in the second, the insistence that all allowed transitions minimised 
a suitably defined local packing fraction led in fact to an intricate coupling between 
the grains. This physically motivated interaction was extremely nonlocal as well 
as directional. In this way we were able to generate a model that, despite being 
one-dimensional, has an extremely complex ground-state structure, depending on 
the regularity of the grain shape. Our present investigations, to be published else- 
where concern the effect of zero-temperature and finite-temperature tapping 
of this system; our preliminary studies indicate that for regularly shaped grains, 
strong metastability in the achievable ground states is observed. For irregularly 
shaped grains, as in reality, a far better packing is achievable, since orientations of 
irregularly shaped grains are much better able to fill space |Q . 

It is however a rather salutary exercise to see that despite the relative sophis- 
tication of the new model in its inclusion of non-trivial interactions, most of the 
qualitative behaviour of the packing fraction as a function of steady as well as 
annealed tapping, remains unchanged. We expect that quantitative features such 
as two-time correlation functions will be far more non-trivial in the second model 
than the first, although we expect their overall features to be rather similar. It is 
tempting to speculate that the directionality due to gravity (which leads to strongly 
non-Hamiltonian behaviour, since grain couplings are propagating down the pile) 
which unites both first and second models might well be the most important ingre- 
dient that is needed to describe such lattice-based models of shaken sand. 
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